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Abstract 
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D ' We present a new proof of the an Ergodic Theorem for Wide-Sense Stationary Random 

^ . Processes added with a new canonical Sampling Theorem for finite time duration signals 

• ^ | in the frequency domain (periodograms) which is free from the Nyquist interval Sampling 

>V 

restriction. We point out the usefulnees of such theorem in the context of a model of 
random vibrations transmission (pressure fluctuations). 

<N : 
> ' 

m : 
o . 

: 1 Introduction 

A basic problem in the applications of stochastic processes is the estimation of a signal x(t) in 
the presence of an additive interference f(t) (noise). The available information (data) is the 
sum S(t) = x(t) + f(t) and the problem is to establish the presence of x(t) or to estimate its 
form. The solution of this problem depends on the state of our prior knowledge concerning 
of the noise statistics. One of the main results on the subject is the idea of maximize the 
output signal-to-noise ratio ([1]), the matched filter system. One of the most important results 
on the subject is that if one knows a priori the signal in the frequency domain X(w) and, 
mostly important, the frequency domain expression for the noise correlation statistics function 
Sff(w), one has, at least in the theoretical grounds, the exactly expression for the optimum 



filter transference function H opt (w) = k(X* (w))(S '/ f(w))~ 1 ■ e _M "*, here t denotes a certain time 
on the observation interval process [—A, A], supposed to be finite here. 

As a consequence, it is important to have estimators and analytical expressions for the 
correlation function for the noise, specially in the physical situation of finite-time duration 
noise observation. Another very important point to be remarked is that in most of the cases 
of observed noise (as in the turbulence research ([2]); one should consider the noise (at least in 
the context of a first approximation) as an Ergodic Random Process. 

We aim in this note to present in Section 2 [of a more mathematical oriented nature] a rigor- 
ous functional analytic proof of an Ergodic Theorem stating the equality of time-averages and 
Ensenble-averages for wide-sense mean continuous Stationary Random Processes. In Section 3 
somewhat of electrical Engineering oriented nature, we present a new approach for sampling 
analysis of noise, which leads to canonical and invariant analytical expressions for the Ergodic 
noise correlation function Sff(w) already taking into account with the finite time-observation 
parameter which explicitly appears in the structure formulae, a new result on the subject, since 
it does not require the existence of Nyquist critical frequency on the sampling rate, besides of 
removing, in principle, the aliasing problem in the computer-numerical sampling evaluations. 
In Section 4 we present an mathematical-theoretical application of the above theoretical results 
for a model of Turbulent Pressure fluctuations. 

2 A Rigorous Mathematical proof of the Ergodic theo- 
rem for Wide-Sense Stationary Stochastic Process 

Let us start our section by considering a wide-sense mean continuous stationary real-valued 
process {X(t), — oo < t < oo} in a probability space {fl, dfi(X), A G Q}. Here Q is the event 
space and d/j,(X) is the underlying probability measure. 

It is well-know ([1]) that one can always represent the above mentioned wide-sense stationary 
process by means of a unitary group on the Hilbert Space {L 2 (f2), dfi(\)}. Namely [in the 
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quadratic-mean sense in Engineering jargon] 

/+oo 
e iwt d(E(w)X(0)) = e tHt (X(0)) (1) 
-oo 

here we have used the famous spectral Stone-theorem to re-write the associated time-translation 
unitary group in terms of the spectral process dE(w)X(0), where H denotes the infinitesimal 
unitary group operator U(t). We have supposed too that the a- algebra generated by the X(t)- 
process is the whole measure space Q, and X(t) is a separable process. 

Let us, thus, consider the following linear continuous functional on the Hilbert (complete) 
space {L 2 (Q), dfi(X)} - the space of the square integrable random variables on Q 

L(Y(X)) = hm ± dtE{Y{X)X^X)}. (2) 

By a straightforward application of the R.A.G.E. theorem ([3]), namely: 

L(Y(X)) = JjKX)Y(X) I Km ^ dt e~ iwt dE(w)X (X) }=Mmi j\{Ye~ m X}dt 



= [ dfx(X)Y(X)dE(0)X{0,X) 
Jn 



= E{Y(X)P KCI(H) (X(0,X))} (3) 

Here Pkct(h) is the (ortoghonal projection) on the kernel of the unitary-group infinitesimal 
generator H (see eq (1)). 

By a straightforward application of the Riesz-representation theorem for linear function- 



al on Hilbert Spaces, one can see that P Kei ( H )(X (0))dfj,(X) is the searched time-independent 
ergodic-invariant measure associated to the ergodic theorem statement, i.e.: For any square 
integrable time independent random variable Y(X) G L 2 (Q,d(x(X)), we have the ergodic result 
(X(0,\)=X(0). 

tl^fl dtE { X ( t W}= E { P Kcr(H)(X(0))Y} (4) 
J — T 

In general grounds, for any real bounded borelian function it is expected the result (not proved 
here) 

T^oh l dtE if( X ^ Y } = E{P Kc r ( H)f(X(0))-Y} (5) 
J — T 



For the auto-correlation process function, we still have the result for the translated time ( 
fixed (the lag time) as a direct consequence of eq(l) or the process' stationarity property 

t^^t r*^W*W* + 0} = ^W°WO} (6) 

J — T 

It is important remark that we still have the probability average inside the ergodic time- 
averages eq(4)-eq(6). Let us call the reader attention that in order to have the usual Ergodic 
like theorem result - without the probability average E on the left-hand side of the formulae, 
we proceed [as it is usually done in probability text-books ([1])] by analyzing the probability 
convergence of the single sample stochastc-variables below [for instance] 



IT 



T 



Vt = t^ f(X(t))dt (7) 



T 



RAO = 7^ J^dtx(t)x(t + c) (8) 

It is straightforward to show that if E{f(X(t))f(X(t + Q)} is a bounded function of the 
time-lag, or, if the variance below written goes to zero at T — > oo 

4= lim -L I dh [ dt 2 [E{X{t l )X{t 1 + C)X{t 2 )X{t 2 + Q)- 

T^oo 41 1 J_ T J_ T 

- EiX^Xit, + C)}E{X(t 1 )X(t 1 + C)))] = 0, (9) 

one has that the random variables as given by eq(7)-eq(8) converge at T — > oo to the left- 
hand side of eq(4)-eq(6) and producing thus an ergodic theorem on the equality of ensemble- 
probability average of the wide sense stationary process {X(t), —oo < t < oo} and any of its 
single-sample {X(t), — oo < t < +00} time average 

T^L^f f^f^)) = E{PK C r(H)(f(X ))} = E{f(X(t))} (10) 

tToo f T dtY ^^ + = E{X(0)X(C)} 

= E{X(t)X(t + Q} 

= Rxx(C) (11) 
The above written formulae will be analyzed in next Electrical Engeneering oriented section. 



3 A Sampling Theorem for Ergodic Process 

Let us start this section by considering F(w) as an entire complex-variable function such that 

\F(w)\ < ce AH , for w = x + iy G C = R + iR (12) 
and such its restriction to real domain w = x + iO satisfies the square integrable condition 

\F(x)\ 2 dx < oo. (13) 



By the famous Wienner Theorem (ref. [4]) there exists a function f(t) G L 2 (—A, A), 
vanishing outside the interval [—A, A], such that 



1 ^ 



F(w) = I dtf(t)e ltw (14) 



A 

The time-limited function f(t) can be considered as a physically finite-time observed (time- 
series) of an ergodic process on the sample space Q = L 2 (R) as considered in section I (since 
Rxx(0) < +oo). In a number of cases, one should use interpolation formulae in order to 
analyze several statistics aspects of such observed random process. One of the most useful 
result in this direction, however based in the very deep theorem of Carleson on the pontual 
convergence of Fourier series of a square integrable function ([1], [5]) is the well-known Shanon 
(in the frequency-domain) sampling theorem for time-limited signals, namely ([1], [6]) 

v > \ T ) T((w — 

n=— oo v V 1 / 

Note in eq(15) however, the somewhat restrictive sampling condition hypothesis of the Nyquist 
interval condition T > A. Althought quite useful, there are situations where its direct use 
may be cumbersome due to the Nyquist condition on the signal sampling besides the explicitly 
necessity of infinite time observations {^?}„ e z- 

At this point, let us propose the following more invariant sampling - interpolation result of 
ours. Let f(t) be an observed time-finite Random ergodic signal (sample) as mathematically 
considered in eq(12)-eq(14). It is well-known that we have the famous uniform convergent 
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Fourier Expansion for t = A sen 9, with < 9 < +— for the Fourier Kernel in terms of Bessel 

functions 

e+iw t = e+iw (Asene) = J- j n ( w A)e in6 (16) 

n=— oo 

As much as in the usual proof of the Shanon results eq(15), we introduce eq(16) into eq(14) 

and by using the Lebesgue convergence theorem, since f(t) G L 1 ^— A, A) either, we get the 

somewhat (canonical) interpolating formula without any Nyquist - like restrictive sampling 

frequency condition on the periodogran F(w) 

+oo ( r / 2 ] 

F(w)= V AJ n {wA){ / d9 e me f(A sen 9) cos 9 \ (17) 

n=-oo J 

It is worth call attention that the sampling coefficients as given by eq(17) are exactly the 
values of the Fourier transform of the signal g{9) = f(Asen9) cos 9 for — — < 9 < -, i.e.: 



2 2' 



1 r /2 

— / ^e me /(^sen^)cos^ (18) 

2tt y_ 7r/ 



T/2 

1 r A N / 1 / 

= A dt ein arc sin(i) /(*) - E (ai ex p( in arc sen (i)) • /(**) 

which may be easily and straightforwardly evaluated by FFT algorithms, from arbitrary - finite 
on their number - chosen sampling values {/(Asen^)}^^...^ with Oi = arsen (^) [here U 
are the time observation process signal]. Monte-Carlo integration techniques ([7]) can be used 
on approximated evaluations of eq.(18) too. Note that f(w) is completely determined by its 
samples eq(18) taken at arbitrary times U E [—A, A]. 

In the general case of a quadratic mean continuous wide-sense stationary process {X t ,0 < 
t < oo} ([1]), we still have the quadratic mean result analogous to the above exposed result for 
those process with time-finite sampling ([1]) 



/+oo 
e iwt dX w (19) 
■oo 

with the spectral process possesing the canonical form (in the quadratic mean sense) 

X(w)= Jn{wA){ / e m6 d(X t=Ascne )\ (20) 



Let us expose the usefulness of the sampling result eqs(17)-(18) for the very important 
practical engineering problem of estimate the self-correlation function of a time-finite observed 
sampling of a Ergodic process already taking into consideration the time-limitation of the 
sampling observation in the estimate formulae. We have, thus, to evaluate the (formal) time- 
average with time-lag ( > 

i2//(C)= Bm ^ £ T dtf(t)f(t + Q. (21) 

Since we have observed the sampling-continuous function {f(t)} only for the time-limited 
observation interval ([—A, A]), it apears quite convenient at this point to re-write eq(21) in 
the frequency domain (as a somewhat generalized process) (see [1]) where the T — > oo limit is 
already evaluated 

Rff(C) = ^J + °° d Wl dw 2 e^^lirn^ ±, j\e it ^ +w ^F(w 1 )F(w 2 )^ 

dw\F(w)\ 2 e iw<: (22) 



2tt 

It is worth point out that eq(22) has already "built" in its structure, the infinite time-ergodic 
evalution < lim — / (...) >, regardless the time-limited finite duration nature of the observed 



T^oo T 



T 



sample f(t)[9(t + A) - 9(t - A)}. 

Now an expression for eq(22), after substituting eq(17) on eq(22) [without the use of some- 
what artificial aliased A-periodic extensions of the observed sampling f(t), as it is commonly 
used in the literature ([6])], can straightforwardly be suggested possesing the important prop- 
erty of already taking into account (in an explicitly way) the presence of the observation time 
A in the result for the self-correlation function in the frequency domain 

+oo +oo 

'27T 



^//W = -7^E E {[ A2 (9n9 m )} x J n (wA)J m (wA)} (23) 

with 



n=— oo m=— oo 



+oo 

iwC, 



R„{Q = -j=l dwe™*S f ,(w) (24) 



The above written equations are the main results of this Section 3. 
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It is worth call the reader attention that after passing the random signal f(t) by a causal 
linear system with transference function Hff(w) ([6]), one obtain the output self-correlation 
function statistics in the standard formulae 



S w (w) = \H vf {w)\ 2 xS ff (24) 

Finally, let us comment on the evaluation of eq(23) on the time-domain. This step can 
be implemented throught the use of the well-known formula for Bessel functions (see [8] - eq 
6.626), analytically continued in the relevant parameters formulae 

POD 

/ dx e~ ax J[ l (bx)J 1/ (cx)dx = I^,u(a) 
Jo 



r(i/ + i) v ' 1^ ew^ + e + i) 
x + (-0} 

We have, thus, the result: 

1 /" + °° t t 

-= / dw e iwt J n (A W ) ■ J m ( W A) = [(-l) n+m I n , m {--r) + /n,m(-r)] 
V27T J-oo A A 



(25) 



(26) 



where 



In ' m[t) ~ r(m + l) 1 } 



x 



Another point to be called the attention and related to the integral evaluations of Fourier 
Transformed of Bessel Functions are the recurrence set of Fourier Integrals of Bessel functions 
below 

»+oo 



/+oo f poo 

dte iwt J n (w) = (1 + (-l) n ) <^ / dtcos(wt)J n (w) 
oo WO 

+ *(l + (-l) n+1 ) jjf dt sen(wt)J n {w)^ (28) 
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with here, the explicitly expressions T n (t) = 2~ n [(t + Vt 2 - l) n + (t - Vt 2 - l) n ]. 



/•oo 

/ cos(wt) J n (w)dw = < 
Jo 



(-lf—L^ T n (t), for < t < 1 and n = 2k 
0,for < 1 < t 

Jp 00 cos(wt)^ J n _ 1 (w)dw — J °° cos(wt)J n _ 2 (w)dw for n = 2k + 1 



(29) 



sen(wt) J n {w)dw = < 



(_l)*-^_r n (t), for < t < 1 and n = 2A; + 1 
0, for < 1 < t 



2(2k + 1) J °° s-^j^w) - J °° sen(wt) J n _ 2 (w)dw for n = 2k + 2 

(30) 



4 A Model for the Turbulent Pressure Fluctuations (Ran- 
dom Vibrations Transmission) 

One of the most important studies of pressure turbulent fluctuations (random vibrations) is to 
estimate the turbulent pressure component transmission inside fluids ([9], [10]). In this section 
we intend to propose a simple analysis of a linear model of such random pressure vibrations. 

Let us consider a infinite beam backed on the lower side by a space of depth d which is 
filled with a fluid of density p 2 and sound speed v 2 - On the upper side of the beam there 
is a supersonic boundary- layer turbulent pressure P(x,t). The fluid on the upper side of the 
beam which is on the turbulence steadily regime is supposed to have a free stream velocity U^, 
density p\ and sound speed v\. 

The effective equation governing the "outside" pressure in our model is thus given by (d < 
z < oo) 

( Q Q\ 2 ( #2 Q2 \ 

with the boundary condition (2 nd Newton's law) 
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Here the beam's deflection W(x, t) is assumed to be given by the beam small linear deflection 
equation 

d A W(x,t) d 2 W(x,t) n . . . , , 

B J } + m = P(x, t) + ( Pl - p 2 )(x, d, t) (33) 

<7 4 £ O z t 

with B denoting the bending rigidity, m the mass per unit lenght of the beam and p(x, z, t) the 
(supersonic) boundary- layer pressure fluctuation of the exterior medium. 

The searched induced pressure p 2 (x, z, t) in the fluid interior medium (0 < z < d) is governed 

by 

d 2 v ( d 2 d 2 \ 

" 12 (x, z, t) - (v 2 ) 2 — + — p 2 (x, z,t) = (34) 



d 2 t ' ' \d 2 x d 2 z / 

dp 2 (x,d,t) d 2 

-^z— = - p2 Wt WM (35) 
The solution of eq.(31) and eq.(32) is straightforward obtained in the Fourier domain 

~ / 7 \ J v^w + U^k) 2 , 

I v/^A; 2 - (ty + U^ky 



x exp x \Jv\k 2 - (w + Uookyj (z - d) (36) 

where the deflection beam W(k, w) is explicitly given by 

W(k, w) = {P{k, w) — p 2 (k, w, d) 

V L.4 2, Pi( w + Uookfvi \ ( 
x < Bk — mw H = > . 37 

\ ^vfk^-iw + u^-kyj 

At this point, we solve our problem of determining the pressure p 2 (k,w,z) in the interior 
domain eq.(34)-eq.(35) if one knows the pressure p 2 (k,w,d) on the boundary z = d. Let us, 
thus, consider the Taylor's serie in the z-variable (0 < z < d) around z = d, namely 

ft(fc, «,, z) = ftflk, „, d) + ^^1^ -<<) + ■■• + ^^f^ - d)» + • • ■ (38) 

From the boundary condition eq.(35), we have the explicitly expression for the second- 
derivative on the depth z 

dp 2 (k,w, z) 



dz 



z=d 



+p 2 w W(k,w) 



+ p 2 w 2 [P(k, w) - p 2 (k, w, d)} x{Bk*- mw 2 + -y==== > ( 39 ) 
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The second z-derivative of the interior pressure (and the higher ones!) are easily obtained 
recursively from the wave equation (34) (k > 0, k G Z + ) and eq(39) 



^Mk. z, w) ^ = x + (v 2 )V) j ^(t, „) j w (40) 

Let us finally make the connection of this random transmission vibration model with the 
section 3 by calling the reader attention that the general turbulent pressure is always assumed 
to be expressible in a integral form 

/+oo r+oo 
dw / dk e lkx e i{w ~ ku)t G{k) ■ F{w) (41) 
■oo J — CO 

where F(w) is the Fourier Transform of a time-limited finite duration sample function of an 
Ergodic Process ([1]) simulating the stochastic-turbulent nature of the pressure field acting 
on the fluid with the exactly interpolating formulae eq.(17)-eq.(18). Numerical studies of this 
random vibration transmission will be presented elsewhere. 
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